Reconciling J/xjj production at HERA, RHIC, Tevatron, and LHC with NRQCD 

factorization at next-to-leading order 
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We calculate the cross section of inclusive direct J/ip hadroproduction at next-to- leading order 
(NLO) within the factorization formalism of nonrelativistic quantum chromodynamics (NRQCD), 
including the full relativistic corrections due to the intermediate X S^\ 3 5*[ 8 ', and 3 Pj 8 ' color-octet 
states. We perform a combined fit of the color-octet (CO) long-distance matrix elements to the 
transverse-momentum (pr) distributions measured by CDF at the Fermilab Tevatron and HI 
at DESY HERA and demonstrate that they also successfully describe the pr distributions from 
PHENIX at BNL RHIC and CMS at the CERN LHC as well as the photon-proton cm. energy and 
(with worse agreement) the inelasticity distributions from HI. This provides a first rigorous test of 
NRQCD factorization at NLO. In all experiments, the CO processes are shown to be indispensable. 



PACS numbers: 12.38.Bx, 13.60.Le, 13.85.Ni, 14.40.Gx 

The factorization formalism of nonrelativistic QCD 
(NRQCD) 1| provides a rigorous theoretical framework 
for the description of heavy-quarkonium production and 
decay. This implies a separation of process-dependent 
short-distance coefficients, to be calculated perturba- 
tively as expansions in the strong-coupling constant a s , 
from supposedly universal long-distance matrix elements 
(LDMEs), to be extracted from experiment. The relative 
importance of the latter can be estimated by means of 
velocity scaling rules; i.e., the LDMEs are predicted to 
scale with a definite power of the heavy-quark (Q) ve- 
locity v in the limit In this way, the theoretical 
predictions are organized as double expansions in a s and 
v. A crucial feature of this formalism is that it takes into 
account the complete structure of the QQ Fock space, 
which is spanned by the states n = 2S + 1 x}'^ with defi- 
nite spin S, orbital angular momentum L, total angular 
momentum J, and color multiplicity a = 1,8. In partic- 
ular, this formalism predicts the existence of color-octet 
(CO) processes in nature. This means that QQ pairs 
are produced at short distances in CO states and subse- 
quently evolve into physical, color-singlet (CS) quarko- 
nia by the nonperturbative emission of soft gluons. In 
the limit v —> 0, the traditional CS model (CSM) is re- 
covered in the case of S'-wave quarkonia. In the case of 
J/tf> production, the CSM prediction is based just on the 
CS state, while the leading relativistic corrections, 
of relative order 0(v 4 ), are built up by the ^q 8 ', 3 «S{ 8 ', 
and 3 pf ] (J — 0, 1, 2) CO states. 

The greatest success of NRQCD was that it was able 
to explain the J/ip hadroproduction yield at the Fermi- 
lab Tevatron while the CSM prediction lies orders 
of magnitude below the data, even if the latter is eval- 
uated at NLO [3, 0]. Also in the case of J/ip photo- 
production at DESY HERA, the CSM cross section at 
NLO significantly falls short of the data 0, Q . Complete 
NLO calculations in NRQCD were performed for inclu- 
sive J/ip production in two-photon collisions Q, e + e~ 



annihilation Q, and direct photoproduction jf|. As for 
hadroproduction at NLO, the CO contributions due to 
intermediate ^q 8 ' and 3 s[ 8 ' states Q were calculated as 
well as the complete NLO corrections to \J production, 
including the 3 sf ] contribution [9j. 

In order to convincingly establish the CO mechanism 
and the LDME universality, it is an urgent task to com- 
plete the NLO description of J/ip hadroproduction by 
including the full CO contributions at NLO, which is ac- 
tually achieved in this Letter. In fact, because of their 
high precision and their wide coverage and fine binning in 
Pt, the Tevatron data on inclusive J/ip production have 
so far been the major source of information on the CO 
LDMEs [lOJ , and the LHC data to come will be even more 
constraining. Previous NLO analyses of J/ip hadropro- 
duction [H 0] were lacking the 3 P| 8 ^ contributions, for 
which there is no reason to be insignificant. This tech- 
nical bottleneck, which has prevented essential progress 
in the global test of NRQCD factorization for the past 
fifteen years, is overcome here by further improving and 
refining the calculational techniques developed in Ref . [6| . 

Invoking the factorization theorems of the QCD parton 
model and NRQCD [lj, the inclusive J/ip hadroproduc- 
tion cross section is evaluated from 

da(AB -t J/ip + X) = ]T J dxdyf l/A {x)f 3/B {y) (1) 

x (a J W[n])da(ij -> cc[n] + X), 

where fi/A{ x ) are the parton distribution functions of 
hadron A, {O j ^ [n]) are the LDMEs, and da(ij -> 
cc[n] + X) are the partonic cross sections. Working in 
the fixed-flavor-number scheme, i and j run over the 
gluon g and the light quarks q — u,d,s and anti-quarks 
q. The counterpart of Eq. ([TJ for direct photoproduction 
emerges by replacing fi/A(%) by the photon flux function 
fy/ e (x) and fixing i = 7. 

We checked analytically that all appearing singularities 
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(4.50 ±0.72) x 10" 2 GeV 3 
(0 J/ ^( 3 S[ 8] )> (3.12 ± 0.93) x 10" 3 GeV 3 
(g) J /^( 3 pW))|(-1.21 ±0.35) x IP" 2 GeV 5 

TABLE I: NLO fit results for the J/tp CO LDMEs. 



cancel. As for the ultraviolet singularities, we renormal- 
ize the charm-quark mass and the wave functions of the 
external particles according to the on-shell scheme, and 
the strong-coupling constant according to the modified 
minimal-subtraction scheme. The infrared (IR) singu- 
larities are canceled similarly as described in Ref. [f|. In 
particular, the 3 Pj 8 ' short-distance cross sections produce 
two new classes of soft singularities, named soft #2 and 
soft #3 terms, on top of the soft #1 terms familiar from 
the S'-wave channels. The soft #2 terms do not factorize 
to LO cross sections; they cancel against the IR singulari- 
ties of the virtual corrections left over upon the usual can- 
cellation against the soft #1 terms. The soft #3 terms 
cancel against the IR singularities from the radiative cor- 
rections to the (C J /^( 3 S'[ 11 )} and (0 J ^( 3 sf ] )) LDMEs. 

We now describe our theoretical input and the kine- 
matic conditions for our numerical analysis. We set 
m c = 1.5 GeV, adopt the values of m e , a, and the 
branching ratios B(J/ip — > e + e~) and B(J/ip — > p + /i~) 
from Ref. and use the one-loop (two- loop) for- 

mula for ai nf \p), with nf = 4 active quark flavors, 
at LO (NLO). As for the proton parton distribution 
functions, we use set CTEQ6L1 (CTEQ6M) [12] at LO 
(NLO), which comes with an asymptotic scale parame- 
ter of Aq^ d = 215 MeV (326 MeV). We evaluate the 
photon flux function using Eq. (5) of Ref. [TIJ with 
the cut-off Qmax = 2.5 GeV 2 15] on the photon vir- 
tually. As for the CS LDME, we adopt the value 
(O j ^( 3 S 1 ^ ] )) = 1.32 GeV 3 from Ref. [ijj. Our de- 
fault choices for the renormalization, factorization, and 
NRQCD scales are p r = /i/ = mx and = m c , respec- 
tively, where mx — + 4m 2 . is the J/tfj transverse 
mass. 

Our strategy for testing NRQCD factorization in J/tp 
production at NLO is as follows. We first perform a 
common fit of the CO LDMEs to the px distributions 
measured by CDF in hadroproduction at Tevatron Run 
II [l7[ and by HI in photoproduction at HERA1 [lj| and 
HERA2 [H[ (see Table Hand Fig.[Q). We then compare 
the p T distributions measured by PHENIX at RHIC [3 
and CMS at the LHC [l]| as well as the W and z dis- 
tributions measured by HI at HERA1 0] and HERA2 
[lj| with our respective NLO predictions based on these 
CO LDMEs (see Fig.©. 

The pt distribution of J / ip hadroproduction measured 
experimentally flattens at px < 3 GeV due to nonper- 
turbative effects, a feature that cannot be faithfully de- 
scribed by fixed-order perturbation theory. We, there- 
fore, exclude the CDF data points with pt < 3 GeV 



from our fit. We checked that our fit results depend 
only feebly on the precise location of this cut-off. We 
also verified that exclusion of the HI data points with 
Pt < 2.5 GeV, which might require power corrections 
neglected here, is inconsequential for our fit. The fit re- 
sults for the CO LDMEs corresponding to our default 
NLO NRQCD predictions are collected in Table Q] In 
Figs. QJa) and (d), the latter (solid lines) are compared 
with the CDF [l]| and HI [Ullll data, respectively. For 
comparison, also the default predictions at LO (dashed 
lines) as well as those of the CSM at NLO (dot-dashed 
lines) and LO (dotted lines) are shown. In order to visual- 
ize the size of the NLO corrections to the hard-scattering 
cross sections, the LO predictions are evaluated with the 
same LDMEs. The yellow and blue (shaded) bands in- 
dicate the theoretical errors on the NLO NRQCD and 
CSM predictions, respectively, due to the lack of knowl- 
edge of corrections beyond NLO, which are estimated 
by varying pL r , pf, and fiA by a factor 2 up and down 
relative to their default values. The /i r , pf, and /i\ de- 
pendencies of a s , the parton distribution functions, and 
the LDMEs, respectively, induced by the renormaliza- 
tion group are canceled only partially, namely through 
the order of the calculation, by linearly logarithmic terms 
appearing in the NLO corrections. Data-over-theory rep- 
resentations of Figs. [IJa) and (d) are given in Figs. (Hb) 
and (e), respectively. In Figs. [IJc) and (f), the default 
NLO NRQCD predictions of Figs. Ufa) and (d), respec- 
tively, are decomposed into their 3 S^\ ^Sq 8 ', 3 sf\ and 
3 Pj 8 ' components. We observe from Fig. QJc) that the 

3 Pj 8 ' short-distance cross section of hadroproduction (ex- 
cluding the negative LDME) receives sizable NLO correc- 
tions that even turn it negative at px > 7 GeV. This is, 
however, not problematic because a particular CO con- 
tribution represents an unphysical quantity depending on 
the choices of renormalization scheme and scale /ia and 
is entitled to become negative as long as the full cross 
section remains positive. Such features are familiar, e.g., 
from inclusive heavy-hadron production at NLO [20| . In 
contrast to the situation at LO, the line shapes of the 
^g 8 ' and 3 Pj 8 ' contributions significantly differ at NLO, 

so that (O j ^ ( x 4 8] )} and (O J ^ ( 3 P(j 8] )) may now be fitted 
independently (see Table HJ. Besides that, the injection 
of HERA data into the fit also supports the independent 
determination of (C J /^( 1 4 81 )} and (O j /^( 3 P { 5 8] )}. Notice 
that (O j /^( 3 Pq 8 ')) comes out negative, which is not prob- 
lematic by the same token as above. In compliance with 
the velocity scaling rules of NRQCD Q, the CO LDMEs 
in Table [J are approximately of order 0(v ) relative to 
(O j /^( 3 S [ i ] )). We read off from Fig. [fla) and (d) that 
the NLO correction (K) factors for hadro- and photopro- 
duction range from 1.30 to 2.28 and from 0.54 to 1.80, 
respectively, in the px intervals considered. 

We observe from Fig. [5] that our NLO NRQCD predic- 
tions nicely describe the px distributions from PHENIX 
(a) and CMS [H| (b) as well as the W distribu- 
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FIG. 1: NLO NRQCD predictions of J/ij> hadro- and photoproduction resulting from the fit compared to the CDF pjj and 
HI [lil. fl5l| input data. The coding of the lines in part (f) of the figure is the same as in part (c). The seeming singularity of 
the 3 fj 8 ' contribution in part (c) is an artifact of the logarithmic scale on the vertical axis.. 



tions from HI flil . Il5j | (c), with most of the data points 
falling inside the yellow (shaded) error band. In all these 
cases, inclusion of the NLO corrections tends to improve 
the agreement. The NLO NRQCD prediction of the z 
distribution (d) agrees with the HI data in the inter- 
mediate z range, but its slope appears to be somewhat 
too steep at first sight. However, the contribution due 
to resolved photoproduction, which is not yet included 
here, is expected to fill the gap in the low-z range, pre- 
cisely where it is peaked; the overshoot of the NRQCD 
prediction in the upper endpoint region, which actually 
turns into a breakdown at z = 1, is an artifact of the 
fixed-order treatment and may be eliminated by invok- 
ing soft collinear effective theory [2l| . We conclude from 
Figs. [T]and [2] that all experimental data sets considered 
here significantly overshoot the NLO CSM predictions, 
by many experimental standard deviations. Specifically, 
the excess amounts to 1-2 orders of magnitude in the 
case of hadroproduction [see Fig. [TJb)] and typically a 
factor of 3 in the case of photoproduction [see Fig. [IJe)]. 
On the other hand, these data nicely agree with the NLO 
NRQCD predictions, apart from well-understood devia- 
tions in the case of the z distribution of photoproduction 
[see Fig. [2][d)]. This constitutes the most rigorous evi- 
dence for the existence of CO processes in nature and the 
LDME universality since the introduction of the NRQCD 



factorization formalism 15 years ago [l|. 

We should remark that our theoretical predictions refer 
to direct J / ip production, while the CDF and CMS data 
include prompt events and the HI and PHENIX data 
even non-prompt ones, but the resulting error is small 
against our theoretical uncertainties and has no effect 
on our conclusions. In fact, the fraction of J ftp events 
originating from the feed-down of heavier charmonia only 
amounts to about 30% [22| for hadroproduction and 15% 
[l5| for photopoprduction, and the fraction of J/ tp events 
from B decays is negligible at HERA [l5[ and RHIC en- 
ergies. 

We thank G. Kramer for useful discussions and B. Ja- 
cak and C. Luiz da Silva for help with the comparison 
to the PHENIX data [l~8j]. This work was supported in 
part by BMBF Grant No. 05H09GUE, DFG Grant No. 
KN 365/6-1, and HGF Grant No. HA 101. 

Note added. At the final stage of preparing this 
manuscript, after our results were presented at an in- 
ternational conference [23| , a preprint [24[ appeared that 
also reports on a NLO calculation of J/tp hadroproduc- 
tion in full NRQCD. Adopting their inputs, we find agree- 
ment with their results for the 3 sf ] , 1 sf\ 3 Sf \ and 3 P l f ] 
contributions. 
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FIG. 2: NLO NRQCD predictions of J/iji hadro- and photoproduction resulting from the fit compared to PHENIX [18], CMS 
and HI [13. Ha] data not included in the fit. 
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